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We present measurements of the differential cross section and A recoil polarization for the 
7P — > K + A reaction made using the CLAS detector at Jefferson Lab. These measurements cover 
the center-of-mass energy range from 1.62 to 2.84 GeV and a wide range of center-of-mass K + 
production angles. Independent analyses were performed using the K + p-K~ and K + p (missing 7r~) 
final-state topologies; results from these analyses were found to exhibit good agreement. These 
differential cross section measurements show excellent agreement with previous CLAS and LEPS 
results and offer increased precision and a 300 MeV increase in energy coverage. The recoil polar- 
ization data agree well with previous results and offer a large increase in precision and a 500 MeV 
extension in energy range. The increased center-of-mass energy range that these data represent will 
allow for independent study of non-resonant K + A photoproduction mechanisms at all production 
angles. 

PACS numbers: 11.80.Cr,11.80.Et,13.30.Eg,14.20.Gk,14.20.Jn,14.40.Aq,25.20.Lj,25.75.Dw 



I. INTRODUCTION 

The 7p — > K + A reaction is a promising channel for 
the study of excited nucleon resonances. Because of the 
pseudoscalar nature of the K + and the self-analyzing de- 
cay of the A baryon, measurement of all polarization 
observables for this channel is experimentally possible. 
Precise measurements of these polarization observables, 
in addition to the unpolarized differential cross section 
(da /d cos Ox™", where O'jf"" is the K + polar angle in the 
center-of-mass frame), will lead to a full characterization 
of the channel and an exciting opportunity to assess the 
contributions of resonant and non-resonant photoproduc- 
tion mechanisms. The channel is further simplified by the 
isospin structure of the final state, which allows coupling 
only to I — \ N* intermediate states and not the I = § 
A* states. 

Previous large-acceptance measurements of the 
7j> — ► K + A differential cross section have been made 
by the SAPHIR [TJ |2J |2j and CLAS 0] collaborations. 
The most recent SAPHIR results [3] are formed from 
roughly 5.2 x 10 4 events and span the center-of-mass 
energy (y/s) range from threshold (1.61 GeV) to 
w 2.4 GeV. The previous CLAS results draw from 
approximately 5.6 x 10 5 events and represent the *Js 
range from threshold to 2.5 GeV. Though differing 
by an order of magnitude in statistics, these results 
do exhibit troubling discrepancies. The SAPHIR re- 
sults are systematically lower (w 20%) than those of 
CLAS at forward angles. While both results show an 
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enhancement in da/d cos O^™- at \fs m 1.9 GeV, this 
enhancement is much more pronounced in the CLAS 
results (especially for cos6*^ m - > 0). Other cross section 
measurements from LEPS at forward |5J and backward 
[5] angles appear to agree with the CLAS results, but 
do not overlap with the regions of the CLAS/SAPHIR 
discrepancy. 

These differences have led to difficulties in interpreta- 
tion of the N* contributions to K + A production. Several 
studies have found evidence for contributions of different 
resonances dependent upon which results are considered. 
The da/d cos#|^ m ' shape discrepancy at s/s ~ 1.9 GeV 
is especially problematic, and partial- wave analyses have 
produced varied explanations for resonant contributions 
in this region, including D i3 [7J, Pi 3 [5], Pn [9], and 
.V M HO] states. 

In this paper, we present measurements of the 
7J3 — > K + A differential cross section and A recoil polar- 
ization (Pa) taken from the CLAS glla dataset. We have 
produced separate analyses using the K + pir~ and K + p 
(missing 7r~) topologies and found these results to be in 
agreement. The K + A measurements presented are the 
most precise to date, and represent an extension of the 
observed yfs range for da/d cos O^" 1 - and P\ of 300 MeV 
and 500 MeV, respectively. These da/d cos #^ m - results 
show agreement with previous CLAS and LEPS results 
and the P\ results agree well with previous world data. 

With several theory groups already pursuing single- 
and coupled-channel partial-wave analyses including the 
7P — > K + A reaction, the results presented herein will 
offer new constraints to pre-existing models. The fine 
center-of-mass-energy binning of these results, especially 
the A recoil polarization, are especially interesting as 
they show previously unseen structure. These results also 
present the first large acceptance measurements of the re- 
action at center-of-mass energies between 2.53 GeV and 
2.84 GeV, an energy regime in which production appears 
to be dominated by non-resonant processes. Previous 
partial-wave analyses have produced non-resonant mod- 
els by constraining only to forward production angle data 
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or by fitting both resonant and non-resonant components 
simultaneously. These da/d cos 9^ m - and Pa data could 
allow for independent study of non-resonant production 
mechanisms at all production angles. 



II. EXPERIMENTAL SETUP 

Data were collected with the CEBAF Large Accep- 
tance Spectrometer (CLAS) located in Experimental Hall 
B at the Thomas Jefferson National Accelerator Facility 
in Newport News, Virginia. The present results are from 
the analysis of the CLAS glla dataset, collected during 
the period of May 17 - July 29, 2004. Photons were pro- 
duced via the bremsstrahlung process using a 4.023 GeV 
electron beam incident on a gold foil. The Hall B tag- 
ger assembly facilitated measurement of the energies of 
recoil electrons using a dipole magnetic field and scintil- 
lator hodoscope; these electron energies were then used 
to calculate the energy of associated photons [IT]. After 
collimation, these photons were incident on the physics 
target, a cylindrical kapton chamber 40 cm in length and 
4 cm in diameter, filled with liquid hydrogen. Measure- 
ments of the target temperature and pressure allowed for 
calculation of the target density with a relative uncer- 
tainty of 0.2%. 

The CLAS detector is composed of tracking and tim- 
ing detector subsystems arranged with six-fold symme- 
try about the beamline (i.e. in six sectors). Trajecto- 
ries of charged particles were deflected by a non-uniform 
toroidal magnetic field with a maximum magnitude of 
1.8 T. Placement of the physics target allowed for recon- 
struction of charged tracks leaving the target at polar 
angles between 8° and 140°. Charged particle tracking 
was accomplished with three sets of wire drift chambers 
per sector. Event timing information was supplied by 
the start counter, a thin, segmented scintillation detec- 
tor placed between the physics target and the innermost 
tracking components, and the time-of- flight (TOF) wall, 
a bank of 48 scintillator bars located beyond the outer- 
most tracking component in each sector. The detector 
subsystems combined to produce an average relative mo- 
mentum resolution of approximately 0.5%. A more de- 
tailed description of the CLAS detector can be found in 
Rcf. [JJ. 

Event triggering required coincident signals from the 
photon tagger and the CLAS level 1 trigger. The CLAS 
level 1 trigger required that two different sectors observe 
a coincidence between timing signals from the TOF and 
start counter scintillators. The signal from the tagger 
consisted of an OR combination of roughly two-thirds 
of the tagger's timing scintillators, which corresponded 
to photons of energy greater than 1.58 GeV. The timing 
scintillators corresponding to lower-energy photons were 
omitted from the trigger in order to reduce the number of 
recorded events generated by photons below the produc- 
tion threshold for many hadronic final states. While the 
number of such events was greatly reduced by this trig- 



ger, events generated by photons with energies between 
wl.O GeV and 1.58 GeV could be recorded due to an ac- 
cidental coincidence with a recoil electron in one of the 
valid tagger elements. For the photon spectrum below 
this energy, a flux renormalization was applied based on 
the probability of such events. With this trigger, physics 
events were recorded and written to disk at a rate of 5 
kHz, only a small fraction of which were relevant to this 
analysis. 



III. DATA AND EVENT SELECTION 

The loose electronics trigger described in the last sec- 
tion allowed for a large number of events to be recorded 
(rj 20 x 10 9 ). Because only a small fraction of these 
events were 7p — > K + A signal, a series of data selec- 
tion cuts was developed to omit events irrelevant to this 
analysis (background). 

Before physics analysis, the dataset was calibrated. 
Timing spectra of the photon tagger, start counter, and 
TOF subsystems were investigated and corrected. Drift 
times from each of the tracking chambers and pulses from 
TOF scintillators were compared and calibrated. After 
these corrections were made, tracks were "reconstructed" 
from raw tracking signals and matched with hits in the 
start counter and TOF detectors. Energy and momen- 
tum corrections were then applied to individual tracks to 
account for imperfections in the magnetic field map and 
detector alignment, and energy losses for particles that 
traveled through the target, detector material, and air. 
Small corrections were also applied to incident photon 
energies to account for slight deformations in the tagger 
hodoscope geometry. 

As the CLAS detector is optimized for detection of 
charged particles, only the charged decay mode of the A 
(A — > pn~) was considered in this analysis. Two separate 
analyses of this reaction were performed: a three-track 
analysis requiring detection of all three of the final-state 
particles, and a two-track analysis requiring only the re- 
construction of K + and p tracks. Possible three-track 
data events were skimmed from the dataset using a 4- 
constraint kinematic fit to the jp — > K + pir~ hypothe- 
sis. Both permutations for the positive track mass hy- 
potheses were tested in all kinematic fits. This fit im- 
posed energy and momentum conservation by varying 
the three-momenta of the detected particles within their 
measurement uncertainties assuming that no undetected 
particles were involved in the event (missing energy and 
three-momentum were constrained to zero for a total of 
four constraints). A probability that each event came 
from the desired reaction (confidence level) was then cal- 
culated from the variations in momenta and the mea- 
surement uncertainties. For the three-track data, events 
with confidence levels less than 1% were removed from 
the analysis. 

Possible two-track data events were selected by per- 
forming a 1-constraint kinematic fit to the "fp — > K + p 
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(missing 7r~) hypothesis and removing events with con- 
fidence levels < 5%. Because the tt~ was not recon- 
structed, this fit imposed only a single constraint that 
the missing mass be that of a tt~ . In order to produce 
results for which uncertainties are dominated by system- 
atic rather than statistical uncertainties, it was sufficient 
to analyze only 28% of the full dataset in producing the 
two-track sample. Both data samples were then sepa- 
rated into 10-MeV-wide y/s bins. The uncertainty in the 
resulting differential cross section measurements due to 
differences in signal lost to the confidence level cuts in 
data and Monte Carlo was estimated to be 3% [IB] . 



IV. BACKGROUND REDUCTION 

Different methods for background subtraction were de- 
veloped for the two- and three-track analyses. The skim 
described in the previous section used a kinematic fit, 
considering tracks' four-momenta and detector resolution 
to select event candidates for the 7p — ► pK + ir~ reac- 
tion. Particle identification was then refined by consider- 
ing timing information for each of the positively charged 
final-state tracks. 

For the three-track analysis, the post-kinematic-fit 
data sample was relatively free of background; back- 
ground events comprised less than 2.5% of the sample 
for all values of center-of-mass energy, \/s. To further in- 
vestigate the nature of this sample, the calculated mass, 
m c , was constructed for each track according to 

m c = v^FW?< (!) 

where p and (3 are the momentum and velocity for the 
particle as calculated from tracking and timing informa- 
tion. By considering two-dimensional histograms of the 
calculated masses of the hypothesized proton and K + 
tracks, the nature of the remaining background is dis- 
cernible (see Fig. [T|. Region (i) in Fig. [TJcontains events 
for which the tracks have calculated masses appropriate 
of the p and K + , thus identifying it as a signal-rich re- 
gion. Region (m) contains events for which the K + track 
is actually a misidentified proton and the proton track is a 
misidentified ir + (i.e., the event is a misidentified ptt + tt~ 
final state) . The majority of the background events lie in 
this region. Region (ii) represents events that passed the 
kinematic fit with the proton and K + tracks reversed 
(i.e., K + misidentified as a proton, proton misidenti- 
fied as a K + ). Events that populate this region are also 
present in region (i) with the correct identification of p 
and K + . Region (iv) contains events for which the pro- 
ton track does not appear to have an appropriate calcu- 
lated mass; however, further investigation of these events 
reveals that they are "fp — > K + pir~ events for which the 
proton timing information was distorted by the detector. 
(This effect is also present in the Monte Carlo, so this 
small fraction of signal events was not removed from the 
analysis.) To remove events from regions (ii) and (m) 



from the analysis, a loose two-dimensional cut on the cal- 
culated masses for the K + and proton tracks was used, 
requiring m c (p) > 0.800 GeV OR m c (K+) < 0.800 GeV. 

In each y/s bin, we fit a Gaussian function to the miss- 
ing mass off of the K + distribution and removed any 
events for which this quantity was greater than 2.5 stan- 
dard deviations from the mean. Finally, fiducial cuts 
were applied to remove events from kinematic or detector 
regions that could not be reliably modeled. The three- 
track data sample included ss 1 .5 x 10 6 signal events occu- 
pying the y/s range from 1.63 to 2.84 GeV, which analysis 
and fiducial cuts reduced to w 6.5 x 10 5 events with less 
than 1% background content at all y/s. Signal loss to 
particle identification cuts for this topology was found to 
be less than 0.11%. 

Because of the less-restrictive kinematic fit for the two- 
track analysis, this data sample had a larger percentage 
of background events. To mitigate this, we first applied 
the same cut on the calculated proton and K + masses 
described above. For y/s < 1.660 GeV, an additional 
two-dimensional calculated mass cut was used to remove 
pn + background. This cut kept events for which 

m c (p) < m c {K + ) + 0.75 GeV/c 2 . (2) 

These cuts remove roughly half of the background events, 
and the Feldman-Cousins method [H] was used to esti- 
mate signal loss to be less than 0.45% for y/s < 1.66 GeV 
and less than 3.4% for y/s > 1.66 GeV. Fiducial cuts were 
applied as in the three-track analysis. 

We then applied to the two-track sample an event- 
based background subtraction technique described in 
Ref . [TS] . This procedure assigns to each event a quality 
factor (Q-factor) that was used to weight the event's con- 
tribution to the fit and the differential cross section calcu- 
lation. We defined a metric based upon the cosine of the 
K + production angle in the cm. frame (cos6*^ m ) and 
cosine of the proton momentum polar angle (cos 0p HF ) 
and proton azimuthal angle (<pp HF ) in the A helicity 
frame. For a given event i, this metric was then used to 
identify event i's 100 "nearest neighbors." An unbinned 
maximum-likelihood fit of a Gaussian signal (s(m)) and 
linear background (6(m)) functions was then performed 
to the missing mass off K + values for these 101 events. 
For event i, the Q-factor was then calculated from the 
signal and background functions: 

Q i = s i /{s i + b i ). (3) 

An example of the signal and background separation in 
a single y/s bin is shown in Fig. [2] We then summed the 
Q-factors for all events to estimate the number of signal 
events present in the two-track data sample after all cuts 
to be w 1.66 x 10 6 . 



V. DETECTOR ACCEPTANCE 

Monte Carlo modeling of the detector acceptance was 
done using GSIM, a GEANT-based simulation of the 
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FIG. 1: (Color On-line) calculated mass of p tracks vs. calculated mass of K + tracks for all events in the three-track data 
sample: Regions (i) and (iv) correspond to jp — ^ K + A signal events, while regions (it) and (Hi) represent background and are 
cut from the analysis. See text for details. 
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FIG. 2: The histograms above show the missing mass off K + distributions for events in the y/s = 2.105 GeV bin of the 
two-track data sample. Figure (a) shows the distribution for events with —0.05 < cos(9^- m ' < 0.05. The unshaded histogram 
shows all events, whereas the shaded histogram shows the same events weighted by (1 — Qi). Figure (b) shows all events in the 
t/s — 2.105 GeV bin (no restriction on cos#^- m ) with the same shading scheme. See text for details. 



CLAS detector. 3 x 10 8 -fp — > K + A events were pseudo- 
randomly generated according to a phase-space distribu- 
tion. GSIM was used to simulate the detector's effects on 
these "raw" events, and a set of "accepted" Monte Carlo 
events was obtained after processing. (GEANT was also 
used to simulate the A — > pn" decay, assuming no net 
polarization for the hyperons.) Corrections accounting 



for the efficiency of the event trigger were applied based 
on efficiencies of individual timing components (TOF and 
start counter). Accepted Monte Carlo events were pro- 
cessed with the same series of analysis cuts as the data 
events. An additional momentum smearing algorithm 
was used to match the momentum resolution of the ac- 
cepted Monte Carlo events with that of the data. More 
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detailed descriptions of the full detector simulation can 
be found in Ref. [TBj . 

In order to form an accurate characterization of 
CLAS's acceptance for the — > K + A reaction, the ac- 
cepted Monte Carlo events were weighted to resemble the 
data following the work in Ref. [13 . To do this, we ex- 
panded the scattering amplitude, M., for the reaction in 
a large set of basis states: 

2 

M (A, a), (4) 

J=\p=± 

where m 7 , rrii, and m\ are the spin projections along the 
beam direction of the incident photon, target proton, and 
A, respectively; X represents the physically significant 
kinematic quantities (cos6*^- m -, cos8p HF and <pp HF ); A 
are the s-channel partial-wave amplitudes for an inter- 
mediate spin-parity J p state (using \ < J < ^ and 
P = ±); and a denotes a set of 34 fit parameters. For 
this expansion, the s-channcl partial-wave amplitudes for 
the 7p — > K + A reaction serve as basis states and were 
calculated for each data and Monte Carlo event using the 
qft++ package [IT]- Estimators, d, for the fit parameters 
were then obtained via unbinned maximum-likelihood fits 
to the data in each y/s bin. We stress that the results 
of this fit are not interpreted as physically meaningful 
(i.e., they do not describe resonant contributions to the 
reaction); the fit results merely express the expansion 
scattering amplitude prescribed by the data. 

Based on this expansion of the data, we then assigned 
to each Monte Carlo event i a weight, Ii, given by 

h= \M m ^ m ^ mA (Xi,a)\ 2 . (5) 

m T ,mi : m\ 

The weighted accepted Monte Carlo matches the data in 
distributions of all physically significant observables and 
their correlations (See Fig. [3]), indicating that our set 
of basis states is large enough to insure a good fit. We 
then calculate the detector acceptance, 77, for a region 
described by kinematic variables X as 

ri{X)=(^l)jl{y j lX (6) 

where the numerator and denominator sums are over the 
accepted and raw Monte Carlo events, respectively. 

It should be noted that this expansion method of ac- 
ceptance calculation replaces the method of attempting 
to generate a physics model for Monte Carlo generation 
via iteration, a method that can be complicated if back- 
ground is present in the data. This expansion method 
also allows for separate calculation of the acceptance for 
each of the six sectors of CLAS. To estimate systematic 
uncertainties in this acceptance calculation, acceptance- 
corrected data yields were calculated independently for 
each sector of CLAS in twelve yfs bins. By considering 
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FIG. 3: (Color On-line) Shown above are cos#^ m ' distribu- 
tions for events in the yfs — 2.105 GeV bin. The data and 
unweighted accepted Monte Carlo are shown by black circles 
and blue crosses, respectively. The Monte Carlo distribution 
weighted according to the scattering amplitude expansion is 
shown by red triangles; this weighted distribution matches 
that of the data events. See text for details. 



the variation in these acceptance-corrected yields in each 
bin, a -^/s-dependent uncertainty (a^yfs)) in the accep- 
tance calculation was determined to be 

0^(71) = 0.0243^/GcV - 0.00890, (7) 

which ranges from 3.0% at threshhold to 6.0% at 
yfs = 2.835 GeV. 



VI. EXTRACTING P A 

The expansion of the data described in Eq. Q allows 
for an elegant and efficient extraction of the A recoil po- 
larization (similar to that described in Ref. [12] )■ In this 
expansion, we have chosen to represent the photon mo- 
mentum as the z direction, choosing the remaining axes 
such that the transverse K + momentum is parallel to x 
and y extends perpendicular to the reaction plane. The 
A recoil polarization, Pa, is a measure of the A baryon's 
polarization out of the reaction plane, i.e. along the y- 
axis. As we have written our amplitudes in terms of the 
z-projections of the photon, target proton, and A spins, 
the recoil polarization at a given value of the kinematic 
variables, X, can be easily projected from the scattering 
amplitude. 

To do so, we first construct a two-component wave 
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function, -0, given by 



MX) = I ^ 



(8) 



where M — ± indicates the spin projection of the A along 
the z-axis and A is the scattering amplitude evaluated for 
the appropriate spin projections and kinematics. Pa is 
then projected with a simple application of a y , the Pauli 
spin matrix in the S'z-basis: 



Pa = Jf E ^ a y^ 



N 



*m T ,m,i ,- 



where 



2V = 



E E 



.M 



P0I S 



(9) 



(10) 



(11) 



is a normalization factor. 

This projection method presents several benefits over 
traditional methods of fitting proton asymmetry distribu- 
tions. Because expansion of the scattering amplitude is 
used for the acceptance calculation method, extraction of 
Pa from amplitudes requires little further analysis. The 
traditional method requires independent fits to proton 
momentum asymmetry distributions in a large number 
of (y/s , cos 6k" 1 ) bins. This compound binning of the 
data can lead to low statistics in some kinematic bins, 
making binned % 2 fits in these bins difficult to interpret 
due to large parameter uncertainties. The method pre- 
sented here is both more efficient, requiring only a single 
global fit in each y/s bin (see Sec. [V]), and more sta- 
ble (with respect to iterations and initial parameter val- 
ues) due to the use of unbinned maximum-likelihood fit- 
ting. Lastly, any Pa measurements given by this method 
are constrained to lie within the physical range, i.e. 
Pa € [1,-1], a feature that is not guaranteed by the 
traditional extraction method. 



VII. NORMALIZATION 

The photon flux during this experiment's run period 
was determined by measuring the rate for electrons inci- 
dent on the photon tagger not corresponding to a trig- 
gered physics event in CLAS. Corrections were made to 
account for the live time of the data acquisition sys- 
tem. Photon attenuation between tagger and physics 
target was studied using a total absorption counter down- 
stream of CLAS. Taking these effects into consideration, 
an energy-dependent total photon flux was calculated 
according to the energy segmentation of the tagger ho- 
doscope. More information on the flux normalization cal- 
culation can be found in Ref. 118. 



Faulty tagger electronics prevented accurate electron 
rate measurement for photons in the energy range 
2.730 GeV < ^fs < 2.750 GeV. Intricacies of the event 
trigger also prevented an accurate flux calculation for the 
y/s bin at 1.955 GeV. Events in this y/s bin could be cat- 
alyzed by photons corresponding to both the triggered 
and un-triggered regions of the tagger. Thus, we present 
no differential cross section results for the y/s =1.955, 
2.735, and 2.745 GeV bins. However, as recoil polariza- 
tion measurements do not depend on the photon flux, we 
do present P\ measurements at these energies. 



VIII. SYSTEMATIC UNCERTAINTIES 

By considering acceptance-corrected yields from in- 
dividual sectors of CLAS, we have estimated a y/s- 
dependent acceptance uncertainty between 3% and 6% 
(see Sec. [V}. Uncertainties due to signal loss to parti- 
cle identification cuts have been estimated to be 0.1% 
for the three-track topology and 0.5% and 3% for the 
two-track topology for bins with y/s > 1.660 GeV and 
y/s < 1.660 GeV, respectively. Uncertainty due to kine- 
matic fit confidence level cuts has been estimated to be 
3% using a study of confidence level distributions for 
the jp —> pTr + 7T~ reaction in this same dataset |13j . 
A 0.5% uncertaintly for the A — > pn~ branching frac- 
tion has been included. Uncertainty in the target length 
and fluctuations in its density contribute 0.2%. Uncer- 
tainty in the flux calculation for this dataset, including 
effects of photon transmission efficiency and live-time 
calculations, has been estimated to be 8% [T3]. Sys- 
tematic uncertainties as they contribute to the two- and 
three-track da/d cos d^" 1 ' measurements are outlined in 
Table |U These individual uncertainties are combined in 
quadrature to yield an overall systematic uncertainty for 
da/dcosOx™ 1 ' measurements of 9%-ll%, dependent on 
topology and center-of-mass energy. 

Because measurement of the A recoil polarization does 
not depend on target characteristics or flux normaliza- 
tion, uncertainties associated with these factors do not 
contribute to the uncertainty in Pa. Our Pa extraction 
method provides no a priori method for calculating the 
associated systematic uncertainty. The effect of accep- 
tance uncertainty on Pa has been studied by considering 
results given by alternate acceptance scenarios [TB] and 
has been estimated to be 0.05 for both two- and three- 
track topologies. 



IX. RESULTS 

A. Differential Cross Section 

For differential cross section (da/d cos 6^-™) and re- 
coil polarization (Pa) measurements, each 10-MeV-wide 
y/s bin was further divided into cos Q'j/ 11 ' bins of width 
0.1. Measurements for these angular bins are reported at 
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Error 


pK + ix 


pK 


-J- / \ 

+ (7T ) 


yfs < 1.66 GeV 


y/s > 1.66 GeV 


Particle ID 


0.1% 


3% 


0.5% 


Confidence Level Cuts 


3% 


3% 


3% 


Acceptance 


3%-6% 


3%-6% 


3%-6% 


Normalization 


8% 


8% 


8% 


Target Characteristics 


0.2% 


0.2% 


0.2% 


A — > pn~ Branching Fraction 


0.5% 


0.5% 


0.5% 


Total 


9%-10.4% 


10%-11% 


9%-10.4% 



TABLE I: Systematic uncertainties in this analysis. The added particle identification cut applied to the two-track analysis at 
low \fs leads to a larger uncertainty. These bins are treated separately. 



the acceptance- weighted bin centroids, the mean of the 
bin range with non-zero acceptance. For each topology, 
in each kinematic (yfs jCOsOjf"") bin, da / d cos djf"" was 
calculated according to 

dcos9 c K mKV ' \T{^s) P d t N A ) 

y(Vs,o c K m -) (12) 

(A cos 6 c K m )r](y/s, cos O^" 1 )' 1 ' 

where A tl pt, and It are the target atomic weight, density, 
and length (respectively), Na is Avogadro's constant, 
T(yfs) is the corrected number of photons incident on 
the target for the given yfs bin, A cos B c ^ a - is the angular 
binning width, and y(y/s, cos O^" 1 ) and rj(y/s, cos 9^ m ) 
are the number of data events and acceptance for the 
given kinematic bin. Differential cross section results for 
both two- and three-track analyses are shown in Figs.[4|6] 
The less restrictive two-track analysis presents measure- 
ments at more kinematic points. In total, the two anal- 
yses present da/d cos B c ^ n - measurements at 2076 unique 
kinematic points. Error bars in these figures represent 
statistical uncertainties from the numbers of data events 
and the Monte Carlo acceptance calculation. 

Several noteworthy features are present in the data. 
For yfs > 1.94 GeV, the forward peak in the data is 
very prominent, and for y/s > 2.4 GeV the forward peak 
dominates the differential cross section, suggesting dom- 
inance of f-channel production mechanisms. In the y/s 
range from 2.4 GeV to 2.65 GeV, we observe a bump in 
the differential cross section at intermediate angles, sug- 
gestive of s-channel production. The scale of this feature 
is small compared to the forward peak, however the fea- 
ture's presence in several ^/s bins is quite interesting as 
production at these energies is considered to be predom- 
inantly t-channel. Above y/s 1.92 GeV, a backwards 
peak is present in the data, and for y/s > 2.39 GeV we 
observe the forward and backward peaks to be separated 
by a relatively flat da/d cos 0j£ m \ This backward-angle, 
high-y^ data presents an exciting opportunity to assess 
it-channel contributions to the reaction. 

Agreement between the two analyses is quantified by 
the relative difference, A, at each kinematic point: 

A(Vi,cos^™-) = , X 4= j ^ f, (13) 



where #2(3) and 02(3) are the result and associated statis- 
tical uncertainty from the two-track (three-track) anal- 
ysis, x is the average of the two results, and a l; (y / s) is 
the acceptance uncertainty. This quantity quantifies the 
difference between the two measurements at a given kine- 
matic point relative to their associated statistical and ac- 
ceptance uncertainties (e.g., A = 1 indicates that the dif- 
ference between two points at a given kinematic is equal 
to the sum in quadrature of their respective absolute un- 
certainties) . We find these relative differences to be nor- 
mally distributed (see Fig. [7]) with mean p = —0.136 
and width a — 0.977, indicating that the two results show 
very little systematic offset and are consistent within sta- 
tistical and acceptance uncertainties. 

This level of agreement between the two analyses leads 
us to produce weighted mean differential cross section 
results according to 

gjcv^cos^™-) = (J2 %)/(E ( 14 ) 

% * j j 

where the sums are over the two analyses and x and 
a represent the measured quantity (here da/d cos 6j^ m -) 
and associated statistical uncertainty. The statistical un- 
certainty on these mean values is then given by 

a 2 (^,cos^" 1 ) = l/o?) ( -2 + \ + — ) , 

(15) 

where the correlation factor, p = 0.28, is due to the 28% 
overlap of the two data samples. For kinematic points 
where only a two-track measurement exists, we use it as 
the mean value and account for the slight offset in the 
two results by scaling its uncertainty by 1 + \p\ = 1.136. 

Comparison of these mean da/d cos 6 c ^ n - results with 
results of previous experiments are worth comment. 
Prior to this analysis, the two highest-statistics studies 
of K + A photoproduction (previous CLAS results [I] and 
SAPHIR 2004 [3]) showed troubling discrepancy. Most 
notably, the previous CLAS differential cross sections 
presented a sizable enhancement at y/s w 1.9 GeV at 
nearly all production angles, whereas the SAPHIR re- 
sults showed a monotonically decreasing da/d cos 6^ m - 
for yfs > 1.75 GeV and cos(9^ m - > -0.15. Though the 
magnitude of the discrepancy between these two analy- 
ses does not exceed 40%, the shape discrepancy has a 
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FIG. 7: (Color On-line) Relative difference of the two- and 
three-track da/d cos 6% m ' results. In red is displayed a fit to a 
Gaussian function yielding the indicated mean (fi) and width 
(a). See text for details. 



large impact on interpretation of K + A photoproduction 
mechanisms. 

Fig. [8] shows the results of this analysis plotted with 
previous high-statistics measurements versus y/s in bins 
of center-of-mass K + production angle. The new CLAS 
results confirm the previous CLAS results at most kine- 
matics, most notably at y/s « 1.9 GeV. These new results 
also show agreement with forward [5] and backward [5] 
measurements from the LEPS experiment, a very illumi- 
nating comparison, as the LEPS results lie at kinematics 
which are typically at the edges of acceptance for the 
CLAS and SAPHIR detectors. We note that these new 
CLAS results are the most precise to date, and extend 
the observed y/s range for this reaction by « 300 MeV. 

The two CLAS results show excellent agreement in 
nearly all of the 120 energy bins, but slight systematic 
discrepancies are present for two specific kinematic re- 
gions. The first region is that of extreme forward K + 
production angles (cosOj/ 71 ' > 0.85). In this region, the 
phase space acceptance extrapolation to kaon angles of 
0° used in the earlier CLAS result was probably less ac- 
curate than the method used in the present analysis. At 
the extreme forward angle, the two measurements are 
only marginally consistent within the respective system- 
atic uncertainty estimates. Also, the CLAS run condi- 
tions for the present dataset had the target offset from the 
center of the detector, thus providing improved forward- 
angle acceptance. 

The second region of discrepancy is the four e nergy 
bins from y/s = 1.715 GeV to y/s = 1.745 GeV. Fig. [9(a) 



shows the very good agreement of the two results in the 
y/s = 1.775 GeV bin (just outside this region). This 
bin is an example of the typically very good agreement 
between the two datasets. In Fig. |9(b)| we present the 
comparison for the y/s = 1.745 GeV bin, that with the 
largest discrepancy of the four bins. These discrepancies 



display a dependence on production angle, beginning at 
cos #^ m - sa 0.2 and continuing to the most forward kaon 
angles. The present results are systematically larger than 
the previous CLAS results at these kinematics, the dif- 
ference between the two being larger than the results' 
quoted statistical and systematic uncertainties. We have 
carefully reviewed both analyses, but have been unable 
to identify problems with either. Thus, we are unable 
to offer unbiased guidance on which data set should be 
preferred for these four energy bins. We can only suggest 
that in this very narrow energy range, the reader exercise 
care when fitting to the CLAS differential cross section 
data. 



B. A Recoil Polarization 



Pa results from the two- and three-track analyses are 
shown in Figs, versus cos dj^ 71 ' in bins of y/s. 

Binning for these results is the same as that used for 
da/d cos Q c ^ n - data. Error bars in these plots represent 
statistical uncertainties. A systematic uncertainty based 
on acceptance uncertainty discussed in Sec. [V]has been 
estimated to be 0.05. In some kinematic areas, differen- 
tial cross section measurements were possible, however 
statistics were too low for a reliable Pa measurement. In 
all, we present measurements at 1708 kinematic points. 

As with the da/dcos9'^ Tn - data, we combine the two- 
and three-track results into a weighted mean result as 
prescribed by Eqs. 



(14 



and (15 1. Fig. 13 shows the 



mean results plotted with previous high-statistics results 
from CLAS [19], SAPHIR [3], and GRAAL 20 . This 
figure shows the new CLAS measurement's increase in 
precision and scope, with a nearly 500 MeV increase in 
y/s coverage at forward angles. The angular resolution 
of this CLAS measurement is unparalleled by any other 
measurement. Comparison between these and existing 
results presents no systematic discrepancies, and several 
structures that are hinted at by previous measurements 
are confirmed by these results. 

Several notable structures are present in the Pa data 
over the y/s range from 1.7 to 2.6 GeV. In the forward 
direction for y/s > 1.9 GeV, where the reaction is known 
to be dominated by i-channel, the recoil polarization is 
relatively featureless with respect to y/s. As one looks 
farther back in production angle, i-channcl mechanisms 
become less dominant and undulations in Pa can be seen. 
As an example, at backward angles, a region of large pos- 
itive A polarization is quite obvious at y/s « 2.0 GeV. As 
one looks forward to intermediate angles, the structure 
remains, but its magnitude is decreased. Several other 
bumps are noticeable in Pa at intermediate angles, in- 
cluding those at y/s » 2.15 GeV and » 2.3 GeV. We 
note that for y/s > 2.1 GeV and very forward angles, 
the recoil polarization remains between —0.5 and —1.0 
indicating a large amount of out-of-production-plane po- 
larization. 
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by open green diamonds, and the LEPS results [5] |B] are shown by open black crosses. 
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FIG. 9: (Color On-line) da / d cos 6>^ m ~ vs. 
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Fig. 9(a) shows results corr espon ding the y/s = 1.775 GeV 
bin of this analysis. Fig. 



9(b) shows those of the y/s 



1.745 GeV bin. Comparisons in this bin are typical of a four- 
bin-wide systematic discrepancy between the two datasets. 
See text for discussion. 



C. Model Comparison 

For first-order interpretation of features in the data, 
we compare the average da/dcos6^ m - and Pa data (as 
prescribed by Eq. 14) to the predictions of severa l co n 

photoproduction 
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temporary models of K + photoproduction. Figs 
show the data and predictions of these models vs. 
bins of cos 6 c j^ n \ 

The Kaon-MAID model [2T] is an isobar model that 
treats non-resonant contributions to the channel as 
i-channel exchanges of K+, K*(892), and #i(1270) 
mesons. Though the Kaon-MAID model is versatile, 
the predictions shown here are from a model fit only to 
SAPHIR data. Resonant contributions to the channel arc 
attributed to the established iV(1650) Su, JV(1710) Pu, 
and iV(1720) P 13 states, as well as a iV(1900) D 13 "miss- 
ing" resonance state necessitated by the enhancement 
of the differential cross section at y/s « 1900 GeV. As 



this model was fit to data of a somewhat limited energy 
range, predictions are only available below y/s = 2200 
MeV. Because it was tuned to the previous SAPHIR 
data, scale agreement between the Kaon-MAID model 
and the present data cannot be expected. However, con- 
clusions can be drawn from comparisons of specific fea- 
tures of the data and the model. 

The second model for comparison is the Regge-Plus- 
Resonance (RPR) model (25] developed by the group at 
the University of Ghent. This model treats non-resonant 
contributions with two Regge-ized i-channel exchanges 
described by a K + Regge trajectory and a K* Regge 
trajectory (both with rotating phases), an elegant de- 
scription requiring only three free fit parameters. As 
Regge models are often considered valid only for small 
exchange momenta, the RPR model was tuned only to 
forward-angle (cos#^ m ' > 0.3) differential cross section 
and polarization data from CLAS and previous high- 
energy data [23]. Resonant contributions in the RPR 
model are the A^(1650) Su, iV(1710) Pu, A^(1720) P l3 , 
and ^(1900) Pi 3 states, as well as a "missing" P/13 state 
with a mass of 1900 MeV. It should not be surprising that 
this model agrees well with the current da/d cos 9'^ m - re- 
sults; agreement between these results and the previous 
CLAS results is satisfactory at most kinematics. 

The final model included here is that of the Bonn- 
Gatchina (BG) group jjj], which is the result of a large- 
scale coupled-channel partial-wave analysis of K + A, 
K + TP, and K°Y, + , pn°, mr + and pr\ photoproduction 
data. It should be noted that the model was constrained 
to 7p — > K + A differential cross section, recoil polariza- 
tion, and beam asymmetry data. This model employs 
the operator expansion method, which projects t- and 
w-channcl amplitudes into s-channel partial waves. Res- 
onant production in the K + A channel is represented by 
significant contributions of the iV(1650) Sn and A^(1730) 
P13 states, as well as two "newly observed" iV(1840) Pu 
and 7V(2170) D 13 states. 

Comparison of these models to the new cross section 
results presents some notable observations. Though the 
Kaon-MAID model displays an almost global scale dis- 
crepancy, it is evident that the model's treatment of 
the cross section at y/s w 1.9 GeV (using a "missing" 
P13 state) is too weak. We also note that the Kaon- 
MAID model overestimates the differential cross section 
for slightly backward angles and y/s > 2.0 GeV. The 
RPR and BG models, as they have been tuned to previ- 
ous CLAS results match the present results well at most 
kinematics. Slight discrepancies exist for the BG model 
at middle angles and y/s w 1.9 GeV and for the BG 
and RPR models at forward angles and y/s « 1.7 GeV 
and y/s > 2.4 GeV. At low y/s, it is possible that 
these discrepancies can be accounted for by re-tuning the 
strengths of s-channcl resonances included. 

One feature of the new cross section results that is not 
reproduced by the models is the slight bump visible at 
cos9 c K m - w 0.0 and y/s~ « 2.1 GeV. The PDG lists several 
N* states with single-star-rated couplings to K + A near 
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this mass, however, a more systematic analysis of the 
data should be performed before associating the feature 
with a given state. 

Agreement of these model predictions and the present 
Pa data is not as good. Recall that previous polar- 
ization data for this reaction were sparse compared to 
the present results. At backward production angles 
(cos6^ m - < -0.15), we see both the Kaon-MAID and 
BG models failing to reproduce the large positive polar- 
ization of the Aat\/ssj 2.0 GeV. At cos6 c K m - w -0.5, 
the models also fail to reproduce the negative A po- 
larization for s/s > 2.2 GeV. At intermediate angles 
(-0.15 < cos6^ m - < 0.35), the BG model reproduces 
the recoil polarization for s/s < 1.85 GeV, however all 
three models fail to reproduce the series of bumps in Pa 
above y/s w 1.85 GeV. As the recoil polarization appears 
to be very sensitive to the nature of the resonances in- 
cluded, as well as interference between resonances and 
between resonances and non-resonant mechanisms, these 
discrepancies could mean that the set of resonances that 
each of these models employs is cither incomplete or in- 
correct. It is worth note that for extreme forward angles, 
only the RPR model seems to accurately describe the 
recoil polarization (though some further tuning of the 
model for 0.6 < cosO^" 1 - < 0.8 is called for), lending 
credence to the Regge-ized meson exchange treatment of 
non-resonant production. 

X. CONCLUSIONS 

In conclusion, these CLAS -fp — » K + A differential 
cross section and A recoil polarization results presented 
here are the most precise to date and offer a significant 
extension of the observed center-of-mass energy range. 
We have presented results from independent analyses of 



the data and found them to demonstrate satisfying agree- 
ment. These analyses provide da/ cos #^ m - and Pa mea- 
surements at 2076 and 1708 kinematic points, respec- 
tively. The da /d cos 6^ m - data show satisfying agreement 
with previous CLAS and LEPS results, while extending 
the observed y/s range by 300 MeV. These results also 
provide overwhelming support for the previous CLAS re- 
sult regarding its discrepancy with SAPHIR results. The 
Pa results presented here agree well with all previous re- 
sults and extend the observed yfs range by 500 MeV. 
These high-precision measurements show a rich structure 
in both observables which present an interesting oppor- 
tunity for interpretation of K + A photoproduction mech- 
anisms. 
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FIG. 10: (Color On-line) Pa vs. cos#^- m ' in bins of y/s. Results from the two-track analysis are represented by closed red 
circles, those of the three-track analysis by open blue triangles. All error bars represent statistical uncertainties only. Horizontal 
and vertical axis scales are common for all plots. Physical limits on Pa are indicated by dashed horizontal lines. 
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FIG. 11: (Color On-line) Pa vs. cos#^- m ' in bins of y/s. Results from the two-track analysis are represented by closed red 
circles, those of the three-track analysis by open blue triangles. All error bars represent statistical uncertainties only. Horizontal 
and vertical axis scales are common for all plots. Physical limits on Pa are indicated by dashed horizontal lines. 
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FIG. 12: (Color On-line) Pa vs. cosO^" 1, in bins of i/s. Results from the two-track analysis are represented by closed red 
circles, those of the three-track analysis by open blue triangles. All error bars represent statistical uncertainties only. Horizontal 
and vertical axis scales are common for all plots. Physical limits on Pa are indicated by dashed horizontal lines. 
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FIG. 13: (Color On-line) Pa vs. ^/s (GeV) in bins of cos#^ m \ Results of this analysis are represented by red circles, previous 
CLAS (McNabb, et al. [T9]) results by blue triangles, SAPHIR 2004 (Glander, et al. [3]) by green triangles, and GRAAL 2007 
(Lleres, et al. [20]) by black squares. Physical limits on Pa are indicated by dashed horizontal lines. 
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FIG. 14: (Color On-line) da/dcosO^" 1 ' (fib) and Pa results us. yfs (GeV) in bins of cos#f/™' plotted with several model 
predictions. Average data points are given by Eq. (14 1. Model predictions are those of Kaon-MAID [5T] (solid green line), the 
Bonn-Gatchina group 9 (dashed blue line), and the RPR model [22] (dashed red line). See text for commentary. 
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FIG. 15: (Color On-line) da / dcos6 c £ n ' (/ib) and Pa results vs. y's (GeV) in bins of cos#f v - m ' plotted with several model 
predictions. Average data points are given by Eq. (14 1. Model predictions are those of Kaon-MAID [5T] (solid green line), the 
Bonn-Gatchina group 9 (dashed blue line), and the RPR model [22] (dashed red line). See text for commentary. 
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FIG. 16: (Color On-line) da/dcos 6 c K m - (/jib) 
predictions. Average data points are given by Eq 

Bonn-Gatchina group [9] (dashed blue line), and the RPR model [22] (dashed red line) 



and Pa results vs. y/s (GeV) in bins of cos 9 c ^ m ' plotted with several model 
Model predictions are those of Kaon-MAID 21 (solid green line), the 

See text for commentary. 
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